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Abstract 

Hadron inclusive spectra in pp collisions are analyzed within the modified quark- 
gluon string model including both the longitudinal and transverse motion of quarks 
Ph! in the proton in the wide region of initial energies. The self-consistent analysis 

I shows that the experimental data on the inclusive spectra of light hadrons like pions 

and kaons at ISR energies can be satisfactorily described at transverse momenta 
not larger than 1-2 GeV/c. We discuss some difficulties to apply this model at 
■ energies above the ISR and suggest to include the distribution of gluons in the proton 

unintegrated over the internal transverse momentum. It leads to an increase in the 
inclusive spectra of hadrons and allows us to extend the satisfactory description of 



^ I the data in the central rapidity region at energies higher than ISR. 

1. Introduction 

m , 

lO ! A rather successful description of various characteristics of hadroproduction pro- 

^ I cesses at not large transfer can be obtained by using the approach for describing the 

^ ■ soft hadron-nucleon, hadron-nucleus and nucleus- nucleus interactions at high ener- 

gies based on the topological 1/N expansion in QCD [1, 2], where is the number 
of fiavours or colours, for example, the quark-gluon string model (QGSM) [3, 4], the 
VENUS model [5], the dual parton model (DPM) [6, 7], the coloured-tube models 
^ I [8, 9] and others. The conventional QGSM and DPM models used the parton distri- 

bution functions (PDF) and the fragmentation functions (FF) integrated over the 
internal transverse momenta of partons. The modification of the QGSM including 
the transverse motion of quarks in the initial hadron has been developed in [10] 
and [11, 12]. It allowed us to describe the inclusive spectra of hadrons produced 
in pp collisions as a function of the Feynman variable x and the hadron transverse 
momentum pt. However, up to now there has not been a self-consistent analysis of 
these spectra within the QGSM in the wide region of initial energies from the ISR 
to the LHC ones. In this paper we present the results of the detailed analysis of 
the inclusive spectra of light hadrons like pions and kaons produced in pp collisions 
at ISR energies within the modified QGSM including the internal transverse mo- 
tion of partons in the initial proton. Then we analyze similar spectra of charged 
hadrons produced in central pp collisions at initial energies from 500 GeV up to 7 
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TeV and compare them with the SppS, Tevatron and latest LHC data. We discuss 
some difficulties to apply the modified QGSM to the description of the latest data. 
To avoid these difficulties we suggest to include the so-called unintegrated gluon 
distributions in the proton by analyzing soft hadron production in pp collisions at 
very high energies. 

As is well known, hard processes involving incoming protons, such as deep- 
inelastic lepton-proton scattering (DIS), are described using the scale-dependent 
PDFs. A distribution hke this is usually calculated as a function of the longitudinal 
momentum fraction x and the square of the four-momentum transfer — —Q^, 
integrated over the parton transverse momentum kf. However, for semi-inclusive 
processes, such as inclusive jet production in DIS, electrowcak boson production 
[13], etc., the parton distributions unintegrated over the transverse momentum kt 
are more appropriate. The theoretical analysis of the unintegrated quark and gluon 
PDFs can be found, for example, in [14, 15]. According to [15], the gluon distribu- 
tion function g{kt) at fixed has the very interesting behaviour at small x < 0.01, 
it increases very fast starting from almost zero values of kt- In other words, g{kt) 
in some sense blows up when kt increases and then it decreases at kt close to 100 
GeV/c. In contrast, the quark distribution q{kt) is almost constant in the whole re- 
gion of kt up to kt ~ 100 GeV/c and much smaller than g{x). This parametrization 
of the PDFs was obtained in [15] within the leading order (LO) and next-to-leading 
order of QCD (NLO) at = 10^ (GeV/c)^ and = 10\GeV/c)^ from the known 
(DGLAP-evolved [16]) parton densities determined from the global data analysis. 
At small values of the nonperturbative effects should be included to get the 
PDFs. The nonperturbative effects can arise from the complex structure of the 
QCD vacuum, see for example [17]- [19]. For example, the instantons are some of 
the well-studied topological fiuctuations of the vacuum gluon fields, see for example 
[20]- [22] and references therein. In particular, it is shown [22] that the inclusion 
of the instantons results in the anomalous chromomagnetic quark-gluon interaction 
(ACQGI) which, for the massive quarks, gives the spin-flip part of it. Within this 
approach the very fast increase of the unintegrated gluon distribution function at 
< A;t < 0.5 GeV/c and = 1 GeV/c is also shown. These results stimulated us to 
assume, that the unintegrated gluon distribution in the proton can be included by 
analyzing also the soft hadron production in pp collisions. We discuss this possibility 
at the end of this paper. 

2. Inclusive spectra of hadrons in pp collisions 
2..1 QGSM 

Let us analyze the hadron production in pp collisions within the QGSM [3, 4] or the 
dual parton model (DPM) [6, 7] including the transverse motion of quarks and di- 
quarks in colliding protons [11, 12]. As is known, the cylinder-type graphs presented 
in Fig.l make the main contribution to this process. The physical meaning of the 
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graph presented in Fig.l is the following. The left diagram of Fig.l, the so-called 
one-cylinder graph, corresponds to the case when two colorless strings are formed 
between the quark/diquark (q/qq) and the diquark/quark {qq/q) in colliding pro- 
tons, then, after their break, qq pairs are created and fragmented to a hadron. The 
right diagram of Fig.l, the so called multi-cylinder graph, corresponds to a creation 
of the same two colourless strings and many strings between sea quarks / antiquarks 
q/q and sea antiquarks/quarks g/g in the colliding protons. The general form for 




Figure 1: The one-cylinder graph (left) and the multicylinder graph (right) for the 
inclusive pp — >■ hX process. 

the invariant inclusive hadron spectrum within the QGSM is the following [3, 4, 23]: 

da _2E* da ^ . . ^ 

where E, p are the energy and three-momentum of the produced hadron h in the 
l.s. of coUiding protons respectively; E*, s are the energy of h and the square of the 
initial energy in the c.m.s of pp; x,pt are the Feynman variable and the transverse 
momentum of h; an is the cross section for production of the n-pomeron chain (or 
2n quark-antiquark strings) decaying into hadrons, calculated within the "eikonal 
approximation" [24], see Appendix; the function (f)n{x,Pt) has he following form [11] 
and [25]- [28]: 

(pn{x,pt)=[ dxi [ dx2'4'n{x,pt;xi,X2) , (2) 

Jx+ Jx- 

where 

^„(x,pt;xi,X2) = F!^^\x+,pt;x^)F^:\x.,pt;x2)/F^:^)iO,pt) + (3) 
+F^:^ (x+, p,; xi)F(;) (x_ , pf, X2)/F(;)) (0, p,) + 
2(n- l)F(;)(x+,Pi;xi)4")(a;_,p,;x2)/F(;))(0,Pi) . 

and x± = 0.5(yaJ^~+~rf ± x),Xt = 2^ {ml + Pt)/s, 

F!r\x±,PuX,^2) = J dXfi^\x±,k,)Gr^, (^yk,;pt)] , (4) 
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Here r means the flavour of the valence (or sea) quark or diquark, f^"'\x',kt) is 
the quark distribution function depending on the longitudinal momentum frac- 
tion x' and the transverse momentum kt in the n-pomeron chain; Gr^hi^, kt')Pt) — 
zDT^h{z, kt'iPt), -Dr- >/i(-z, kt'iPt) is the fragmentation function of a quark (antiquark) 
or diquark of flavour r into a hadron h. We present the quark distribution in a 
proton in the factorized form f^{x,kt) = fT{x)gT-{kt), and choose the kf distribu- 
tion of quarks and the FF in the simple Gaussian form g{kt) = {'jq/n)exp{—jqk'^) , 
Qg^hih) = ilF/'^)exp{-jFkt), where kt^Pt- zkt, z = x±/xi^2- Then, the quark 
functions in the n-pomeron chain will be factorized too [11] 

fi'^\x,kt)=f^-\x)gi-\h) , (6) 

where 

gi^'Kh) = ^exp(-7nA;,^), 7n = ^ ■ (7) 
TT n 

The fragmentation function also reads 

Gr^h{z,kt;Pt) = GT^h{z,Pt)9T^h{kt) ■ (8) 
Then, substituting Eqs.(6) into Eq.(4) we get the following form for Fj'^^: 

F'r\x±.Pt;X^,2) = tf\x^^2)Gr^h{z)in{z,Pt) , (9) 

where the function In{z,pt) reads [12, 11] 

in{z,Pt) = — exp(-7^p^), 7^ = -,7^ 2 ' P^ — ■ (10) 
TT 1 -I- npz^ 7g 

When we take the kt distribution of quarks (diquarks) g{kt) and the FF gq-^h{h) in 
the exponential form 

gikt) = ^exp{-Bqkt) , gq^hih) = -^exp{-BFkt) , (11) 
we get the following form for the function In{z,pt) entering into F^'^\x±, pt; Xi^2) 

F^'^\x^,Pt-X,,2) = fl''\x,,2)Gr^h{z)In{z,Pt) , (12) 

Iizr^\-C f Mbpt)bdb 

^n[Z,Pt) ^nj (^2^2 + 52)3n/2(^2 + 52)3/2 ■ \^'') 

Here Jo{bpt) is the zero order Bessel function and the coefficient C„ is determined 
by the normalization equation 

J £pMz,Pt) = 1 (14) 
4 



Equation (14) is similar to the normalization equation for the function In{z,pt) given 
by Eq.(lO), when the quark distributions and the FF are chosen as the Gaussian 
forms. 

• Gluon distribution in proton . 

As is mentioned in the Introduction, the unintegrated gluon distribution in the 
proton at small values of function of kt, increases very fast when kt increases 

and then slowly decreases, according to the instanton vacuum approach for the 
massive quarks [22] at small ~ 1 GeV/c. The similar behaviour for g{kt) was 
obtained within the NLO QCD calculations at large — lO^GeV/c and = 
lO^GeV/c for the massless quarks [15] and in [29, 30] where the parametrization for 
the gluon distribution at the starting value Ql was found as 

gih, X, Ql) ~ Rlix)k^ eM-Rli^)kt) , (15) 

with the parameters given in [32]. The motivation of the form given by Eq.(15) is 
based on the saturation model [31]-[33]. Unintegated gluon distributions were also 
studied in [14] within the dipole model where the similar qualitative kt dependence 
was shown. On the basis of these results we will try to include the contribution of 
gluons in soft pp collisions. 

• Hadron production in central rapidity region. 

According to the Abramovskiy-Gribov-Kancheli cutting rules (AGK) [34], at mid- 
rapidity only Mueller-Kancheli type diagrams contribute to the inclusive spectrum 
of hadrons. In our approach the function F^^^ is calculated in such a way that in 
the central region {y = 0), when a; ~ and ^ ~ 0, it becomes proportional to n and 
satisfies AK cancellation. Thus, 

oo 

Pq{0,Pt) = (l)q{Q,Pt) ^^n(s) = 9{s / Sq)^ (j)q{0, Pt) , (16) 
n=l 

where ^^(x = 0,pt) depends only on pt, sq — 1 Gev^. Considering the gluons from 
incoming protons, which may split into qq pairs, we get an additional contribution 
to the spectrum 

oo 

Pg{x = 0,pt) = 0s(O,pt) Yin- l)a„(s) = (17) 

n=2 

0ff(O,Pt)(^(s/so)^ -(7„d) , 

where A = 0.12, g=21 mb and and is the cross section of exchange of any number 
cut-pomerons. The quantities 

oo oo 

Y nCTnis) = g{s/so)^ ; J2 '^n{s) = CTj^d (18) 
n=l n=l 

were calculated in [24] within the "quasi-eikonal" approximation [24] . Assuming that 
one of the cut-pomerons is always stretched between valence quarks and diquarks 
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which are not coming from the sphtting of ghions, in Eq. (9) we excluded unity from 
n. Finally, we can present the inclusive spectrum at a; ~ in the following form: 

p(Pt) = Pq(x = 0,pt) + Pg(x = 0,pt) = (19) 

We fix these contributions from the data on the charged particles pt distribution, 
parametrizing them as follows: 

<l)g{0,pt) = Agexp{-bgpt) (20) 
0fl(O,Pt) = Agy/p'texp{-bgpt). 

The parameters are fixed from the fit to the data on the pt distribution of charged 
particles at y = 0: Ag ^ 0.1912 ± 0.0064; bq = 7.24 ± 0.11 (Gev/c)-^ and Ag = 
0.0568 ± 0.002; bg = 3.46 ± 0.02 {GeV/c)-\ 

2.. 2 Hard Scattering 

As will be shown below, the approach suggested above will be able to describe many 
data on the inclusive pt spectra in the central pp collisions at not large values of 
Pt < 2 GeV/c. Therefore, at larger pt wc calculate these spectra, which are due 
to the hard pp collision, within the leading order of perturbative QCD (LOQCD). 
According to the model of hard scattering [45] -[53], the relativistic invariant inclusive 
spectrum of the hard process p + p — )■ h + X can be related to the elastic parton- 
parton subprocess i + j i' + j', where i,j are the partons (quarks and gluons). 
This spectrum can be presented in the following general form [47]-[49]: 

p{x-,Pt)hard = e'^^^^^ = 1] / d^^it I d'^kjt [ dXi [ dxj (21) 
[' dzMx,, ku)fAx„ k,t)-^^^^^Dl^iz)d{s + i+u) 

Jo TT at Z'^ 

where s — (pi +^2)^ — 2pi -pa, P\iP2 are the four-momenta of the colliding protons. 
In the cms. of pp (LHC facility, when p\ = —p^) s — 4£'f, i.e., Ei — -\/i/2; 
t = {Pi - Phf ^ -2pi • Ph = -2{EiEh - piPh); u = (p2 - Phf ^ -2p2 • Ph = 
-2{EiEh -piPh); = 2k(ij)J yG; ^(1,2,^) = ^"^(1,1,^) +P?i,2,/i)' P{i,2,h) are the 

three-momenta of hadrons 1,2, h respectively; are the longitudinal momenta 

(relative to pa) of the partons i or j in the pp c.m.s., z is the fraction of the hadron 
momentum from the parton momentum, fij is the PDF, whereas Dij is the FF, 
s, t, u are the Mandelstam variables for the parton-parton elastic scattering, sec the 
details in [47]. Calculating Eq.(22) we used the PDF, FF and daij{s,t)/dt obtained 
within LO QCD [50]- [53]. 
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3. Results and discussion 



In Figs. (2-5) we illustrate how our approach works for the description of the exper- 
imental data on the inclusive spectra of pions and kaons produced in pp collisions 
when we neglect the gluon distributions in the proton. It allows a satisfactory de- 
scription of these spectra as functions of x and pt, up to pt = 1.4 GeV/ c at the initial 
ISR energies y/s = 23.3 — 53 GeV. The ISR experimental data presented in Figs. (2- 
5) are taken from [35]. The quark (diquark) distribution and the FF as a function 
of the internal transverse momentum kt and kt are chosen in the exponential form 
given by Eq.(ll). We took the same parameters Bq and Bp for all the initial ISR 
energies, namely Bq =4.5 GeV/c for both valence and sea quarks of any flavour and 
Bp =2.8 GeV/c for vr^ and Bp =2.4 GeV/c for mesons. Figures (2-5) show 
that the modified QGSM results in a self-consistent satisfactory description of the 
inclusive spectra of light mesons produced in pp collisions at ISR energies at different 
values of x and pt not larger than 1.5 GeV/c. However, one cannot satisfactorily 
describe similar spectra at the SppS, Tevatron and LHC energies ignoring the gluon 
contribution . 




Figure 2: The inclusive spectrum Eda/d^p [mbGeV'^c'^] of tt^ mesons (left) pro- 
duced in pp collision at ^/s = 23.3 GeV; the similar spectrum but for 7r+ mesons 
(right). 



The result of the fit to data using Eq.(20) on the charged hadron inclusive spectra 
is presented in Figs. (7,8). The long-dash curve corresponds to the quark contribution 
Pq{x = 0,pt) given by Eq.(16), and the short-dash line is the gluon contribution 
Pg{x = 0,pt) (Eq.l7) to the invariant yield d^N/dyd'^pt; the solid curve corresponds 
to the sum of both contributions, see Eqs.(19,20). One can see that the conventional 
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Figure 3: The inclusive spectrum Eda/d^p [mbGeV~^c'^] of tt~ mesons (left) pro- 
duced in pp collision at ^/s = 30.6 GeV; the similar spectrum but for vr"*" mesons 
(right). 
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Figure 4: The inclusive spectrum Eda/d^p [mbGeV^^c'^] of vr^ mesons (left) pro- 
duced in pp colhsion at ^/s = 53 GeV; the similar spectrum but for 7r+ mesons 
(right). 

quark contribution $''(?/ = 0,pt) is able to describe the data up to pt <1 GeV/c, 
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Figure 5: The inclusive spectrum Eda/d^p [mbGeV~^c^] of K~ mesons (left) pro- 
duced in pp collision at ^/s = 44.6 GeV; the similar spectrum but for mesons 
(right). 
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Figure 6: The inclusive spectrum Eda/d^p [mbGeV~^c^] of K~ mesons (left) pro- 
duced in pp collision at ^/s = 53 GeV; the similar spectrum but for mesons 
(right). 



whereas the inclusion of the gluon contribution allows us to extend the range of 
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Figure 7: The inclusive spectrum of the charged hadron as a function of pt (GeV/c) in 
the central rapidity region {y = 0) at ^/s = 540, 900 GeV (top) and ^/s = 1.8, 2.36, 7 
TeV. The data are taken from [36, 37, 38]. 



good description up to 2 GeV/c. At larger values of pt the contribution of hard 
processes is not negligible and one has to take them into account on the basis of 
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Figure 8: The inclusive spectrum of the charged hadron as a function of pt (GeV/c) 
in the central rapidity region (y = 0) at y/s =7 TeV compared with the CMS [38] 
and ATLAS [39] data. 



the perturbative QCD. Figure 7 also shows that, according to the experimental 
data, the shape of the inclusive pt spectrum at y = slowly changes as the initial 
energy increases. Its pt dependence appears harder as the energy increases. It 
is probably due to the jet production which occurs even at not large values of pt 
when ^/s increases. In Fig. 8 the same results as in Fig.7(right, bottom) except the 
solid line are presented. The solid line in Fig. 8 corresponds to the matching of the 
calculations obtained within the soft QCD (SQCD) and perturbative QCD (PQCD). 
The SQCD means our calculations within the QGSM including the gluons in the 
proton, see Eq.(19) In Fig. 9 the ratio of our calculations to the experimental data 
[38] at a/s =900 GeV (top) and y/s =7 TeV are presented. It is shown that the 
discrepancy between theoretical results and the experimental data is less than 30 
percent at pt < 2.GeV/c. 

The effective inclusion of the unintegrated distributions of gluons in the proton 
[29, 30] allows us to extend the satisfactory description of the experimental data 
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Figure 9: The ratio of our calculations of the inclusive pt spectrum as a function of 
Pt (GeV/c) in the central rapidity region [y = 0) to the corresponding experimental 
data at i/i =900 GeV (top) and i/i =7 TeV (bottom). The data are taken from 
[38]. 



on the inclusive spectra of charged hadrons at ?/ = up to ~ 2GeV/c. At 
higher values of the transverse momentum we include hard pp collision calculated 
within LO QCD [45]-[53]. Parton interactions should be considered for describing 
the experimental data. These results are presented in Figs. (7,8) at pt > 2GeV/c. 



4. Conclusion 

Our study has shown that soft QCD or the modified QGSM including both the 
longitudinal and transverse motion of partons in the proton is able to describe 
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rather satisfactorily the experimental data on inclusive spectra of light hadrons like 
pions and kaons produced in pp collisions at the ISR energies and at not large values 
of the transverse momenta pt < 1.4 — l.bGeV/c. These calculations were made in 
the self-consistent way, i.e., the parameters entering into the quark distribution and 
the FF, as a function of the internal transverse momentum, are the same at all the 
initial energies ^/s up to about a few hundred GeV. The conventional QGSM or 
DPM model does not include the distribution of gluons in the proton. However, as 
is well known, at large transverse momenta pt of hadrons the gluons in the proton 
play a very important role in description of the experimental data. Therefore, one 
can assume that the contribution of the gluon distribution in the proton to the 
inclusive spectrum of produced hadrons slowly appears when pt increases and it will 
be sizable at high values of pt- This assumption is also confirmed by the increase 
of the unintegrated gluon distribution in the proton at x ~ as a function of the 
internal transverse momentum kt when kt grows [15, 22, 29]. 

Therefore, to illustrate this hypothesis we fit the experimental data on the in- 
clusive spectra of charged particles produced in the central pp collisions at energies 
larger than the ISR starting from 500 GeV up to 7 TeV by the sum of the quark 
contribution pq given by Eq.(16) and the gluon contribution pg (see Eq.(18)). The 
parameters of this fit do not depend on the initial energy in that energy interval. 
This fit shows that the inclusion of only the quark contribution pq allows us to 
describe the data at y = up to ~ 1 GeV/c, while the inclusion of the gluon 
contribution pg results in a satisfactory description of the data up to pt ~2GeV/c. 
The inclusion of the hard pp collision within LO QCD allows us to describe rather 
satisfactorily the data in the wide region of pt- 
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5. Appendix 

• Quark distributions in proton within QGSM 

Let us present the quark distributions in a proton obtained within the Regge theory 
in Refs.[3, 4, 23]. The distributions for valence u and d quarks in the n chain (see 
Fig.l) read 

ff{x) = C„a;-"«(°)(1 - ^)««(0)-2a.(0)+n-l ^ ^(n)(^) ^ ^(n)(^) ■ {I - x) . (22) 
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There are the following relations for the sea u and d quarks: 

fit (x) - fit (^) = & (^) ; /il (^) - /£ (^) = /i? (^) ■ (23) 

For the charmed quarks cc in a proton we have [12] 

= ^.^"^^(l - 2,)2(a«(0)-a^(0))-a^(0)+n-l (24) 

The distributions of the diquarks ud and uu in a proton read 

/i?(^) - C„,x"«(°)-2"-(°)(l - x)-"«(°)+"-^ ; fi^ix) = fi:\x) • (1 - x) . (25) 

Here aR{0),aN{0) and q;^(0) are the intercepts of the Reggeon, nucleon and ip 
Regge trajectories. As is well known (see, for example, [23]), aR{0) = 0.5 and 
q;7v(0) — —0.5. The value for the intercept of the ip Rcgge trajectory is known not 
so well because Regge trajectories of heavy mesons can be nonlinear as functions of 
the transfer t. For the linear -0 Regge trajectory a^{0) — —2.18, whereas for the 
nonhnear ■0(f) the intercept value can be about zero, a^(0) = [41]. The coefficients 
Cj in Eqs. (22-25) are determined by the normalization condition 

ht\x)dx = 1. (26) 

JO 

• Fragmentation functions of quEirks (diquEirks) to tt mesons within QGSM 

The fragmentation functions of quarks (diquarks) to tt mesons G'q(gg)_^^(2;) = 
have the following forms [23, 40, 42]: 

G'„_,+ (z) = ao(l - ^)"«W+^ ; G^^,-{z) = (1 - z)G^^,+ , (27) 



Gd-^'K+{z) — Gu^T^-{z) , Ga-^T^-{z) — Gu-^T^+{z) , (28) 
Guu-..A^) = ao(l - ^)"«(o)-2aB(o)+A . G^^^,-(z) = (1 - z)G^^^,4z) , (29) 



G^a^^^z) = G^,^,-(z) = ao (l + (1 - z)') (1 - ^)"H(o)-2aB(o)+A (33) 

and 

(z) . (31) 

In the above equations aR{0) — 0.5, q;b(0) = —0.5, oq = 0.65 and A = 2q;^(0) < 
p| >^ 0.5, q;r(0) and < Pt > being the slope of the Regge trajectory and the mean 
value of the Tahoe transverse hadron momentum squared. 
• FF for mesons within QGSM 
The FFs for mesons are [43] 

Gu^K+iz) = Gd^Ko{z) - aK{l - ^)-"^(o)+^(l + a,Kz) , (32) 
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Gu^K-{z) — Gfi^jio{z) — G^^x+{z) — Gu^K"iz) — Gj_^K-{z) — (33) 

Gu-^K-{z) = Gu^K+{z) , Gu^K+{z) = Gu^K-{z) 1 (34) 
Gd--fK-{z) — Gd^K+{z) , G^_^K+{z) — Gd-^K-{z) , 

Gs^K^iz) = Gs^ii+{z) = Gs^K-{z) = G,^Ko{z) = (35) 
^^^1-00(0) j-]^ - 2)-"«(o)+^ + axil - 2)-"fl(o)+^+2(i-«^(o)) 

G,^K-{^) = G,^Ko{z) = G,^kA^) = G,^Ko{z) = a^(l - ^)-"«(o)+A+2(i-a,(o))(36) 

Guu-.kA^) = axil - ^)2"«(0)-a.(0)-2a.(0)+A(i ^ ^ (37) 

Gua-.Mz) = ^(1 - ^)^'^«W-^W-^'^-W+^(l + a^K^ + (1 - ^)^) , (38) 

Guu-.K-{^) = Guu^K0{z) = G,^^K0{z)aK{l - ;,)-.(0)-2a.(0)+A+2 ^ ^39) 

Gud-.K-{Z) - G^^^o{z) = ^(1 - ^)-.(0)-2a.(0)+A+2(;L + (1 - ^)2) , (40) 

G^a-.K<z) = ^(1 - ;,)-.(o)-2«.(o)+A+2(2 + «^^^) , (41) 

Here A = 2a'pt - 0.5, 0^(0) = 0.5, ajv(O) ~ -0.5, aiK = 2, 021^ = 5, bx ^ 0.4. 
• FF for p and p within QGSM 

Now, the FFs for protons p and antiprotons p are [44] 

Gu-,p{z) = ap{l - ^)"H(o)-2a.{o)+A . G^_^^-{z) = (1 - ^)G'„^,+ , (42) 

Ga^p{z) = |(1 - z)--(o)-2--(o)+^(3 - z) ; (43) 

Guu->p{z) = Giuu-^p{z) + G2uu-^p{z) , (44) 

G2..^p(^) = a,{l - ^)-««(0)+A+4(l-a.(0)) ^ (4g) 
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,2(ai{(0)-Qjv(0))/'i _ ^\-a«(0)+A 



(47) 




(48) 



(49) 



Guu-.p{z) = G^a-.p{z) = a,(l - ^)"H(0)-2a.(0)+2(l-a.(0))+A 



(50) 




(51) 



Here ap(0) = 0.9 and ap(0) = 0.07. 
• Hcird scattering 

Present now the scheme for the calculations of the differential cross section of 
the parton-parton elastic scattering within the LO QCD. Figure 9 (top) illustrates 
the parton-parton scattering within the LO QCD. In Fig. 9 (middle and bottom) the 
differential parton-parton cross sections calculated within the LO QCD are presented 
[50]. Here A is the chromodynamic constant which has taken from [52, 53], Uf 
is the number of flavours, Ai is the dimensional coefficient to get the dimension 
mb/(GeV/c)^ for dcjij/dt. The four-momentum transfer squared is related to 
the Mandelstam variables s, t, u for the elastic parton-parton scattering [49] 
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